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■ We analyze the mixing of the scalar glueball with the scalar-isoscalar quarko- 

nia states above 1 GeV in a non-local covariant constituent approach. Sim- 
ilarities and differences to the point-particle Klein-Gordon limit and to the 
quantum mechanical case are elaborated. Model predictions for the two- 
photon decay rates in the covariant mixing scheme are indicated. 
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O ■ The possible existence and observable nature of the lowest-lying scalar glueball is cur- 

Q . rently under intensive debate. The discussion on possible evidence for the emergence of the 
"Th ! glueball ground state in the meson spectrum has dominantly centered on the scalar-isoscalar 
PhI resonance /o(1500) [1]. This state has been clearly established by Crystal Barrel at LEAR [2] 
p ' in proton-antiproton annihilation, and is also seen in central pp collisions [3] and J/\l/ decays 

[4]. The main interest in the /o(1500) as a possible candidate with at least partial glueball 
content rests on several phenomenological and theoretical observations. The /o(1500) is 
produced in gluon rich production mechanisms, whereas no signal is seen in 77 collisions 
^ I [5]. Lattice QCD, in the quenched approximation, predicts the lightest glueball state to be 

I a scalar {J^^ = 0^+) with a mass of 1611 ± 30 ± 160 MeV [6]. Although the decay pattern 
of the /o(1500) into two pseudoscalar mesons might be compatible with a quarkonium state 
in the scalar nonet of flavor structure nn = -^{uu + dd), the observed, rather narrow total 
width of about F ^ 110 MeV is in conflict with quark model expectations of about F > 500 
MeV [1,7]. 

There are mainly two major theoretical schemes, each of them split in a variety of 
sub-scenarios, which try to extract and point at the features of the scalar glueball embed- 
ded in the scalar meson spectrum. The first and original ones [1,8] consider the possi- 
ble nonet of scalar mesons above 1 GeV which are overpopulated according to the naive 
quark- ant iquark picture. Experimentally one identifies a surplus state among the observed 
K*{U30), ao(1450), /o(1370), /o(1500) and /o(1710) resonances. The last three isoscalar fo 
states are considered to result from the mixing of two isoscalar qq states and the lowest- 
lying scalar glueball. A recent phenomenological analysis concerning the hadronic decays 
of /o(1370), /o(1500) and /o(1710) into pairs of pseudoscalar mesons is consistent with the 
minimal three-state mixing scenario [9]. The lower lying resonances /o(980) and ao(980) are 
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not taken into account since they presumably are represented by four-quark configurations 
with a strong couphng to virtual KK pairs [10]. 

In a second scenario [11] /o(980) and ao(980) are interpreted as quark-antiquark states, 
which together with the /o(1500) and K*{im) form a nonet. The /o(980) and /o(1500) 
are the isoscalar mesons of the nonet with large mixing as in the pseudoscalar case. The 
broad resonances /o(400 — 1200) (or a) and /o(1370) are considered as one state which is 
interpreted as the light glueball. A similar approach, which allocates most of the glueball 
strength in the /o(1370), is considered in Ref. [12]. Some works [13] in QCD sum rules also 
suggest that the glueball component is present in the low mass /o state. However, a recent 
compilation [14] of theoretical arguments in comparison with experimental signatures seems 
to favor the first scenario, where the glueball strength is centered in the /o(1500) region. 

The aim of the present work is to follow up on the three-state mixing scenario, where the 
glueball intrudes in the scalar quark-antiquark spectrum at the level of 1.5 GeV giving rise 
to the observed /o(1370), /o(1500) and /o(1710) resonances. Up to now many studies in this 
direction have been performed ( [1,7,9,15] and Refs. therein), considering a mixing matrix 
linear in the masses and applying a SO{?>) rotation to find the orthogonal states, which are 
then identified with the three resonances mentioned above. The mixed and unmixed states 
are linked by an 5*0(3) matrix M from which one can read the composition of the resonances 
in terms of the "bare" unmixed states. In the present work we develop a framework, where 
the glueball-quarkonia mixing is described in a covariant field theoretical context. The 
present analysis proposes a way to define the mixing matrix M which allows a comparison 
with the Klein-Gordon case and with the quantum mechanical limit. The resulting matrix 
M is in general not orthogonal which can be traced to the non-local, covariant nature of the 
bound-state systems. 

Work along these lines has been performed in the context of effective models [16-18], 
where the glueball degree of freedom is introduced by a dilaton [17,18] or by a scalar field 
constraint by the QCD trace anomaly [16]. In the present approach we resort directly to the 
elementary gluon fields in a non-local interaction Lagrangian, where the non-locality allows 
to regularize the model and to take into account the "wave function" of the gluons confined in 
the glueball. We also introduce effective non-local quark-quark interactions, working in the 
SU{3) flavour limit, in order to describe the bare, unmixed quark-antiquark states. Mixing 
of these configurations is introduced on the basis of the fiavor blindness hypothesis and the 
glueball dominance, where mixing is completely driven by the glueball. Another interesting 
feature of our study is that we obtain breaking of fiavour blindness, when considering the 
composite nature of the mesons. 

As an elementary input we have to specify the non-perturbative quark and gluon prop- 
agators, which should not contain poles up to ~ 1.8 GeV, thus avoiding the on-shell man- 
ifestation of quarks and gluons (confinement). Many considerations and characteristics are 
independent on the form of the propagators, as will be shown in the numerical results. 

The paper is organized as follows: in Section II we present the details of the model as 
based on a non-local covariant Lagrangian formalism. Section III is devoted to the implica- 
tions of the model, like mixing properties, the definition of the mixing matrix M and the 
comparison with other mixing schemes. We furthermore specify the two-photon decay rates 
of the mixed /o states. In Section IV we present our results and the related discussions. 
Finally, Section V contains our conclusions. 
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II. THE MODEL 



We will consistently employ a relativistic constituent quark/gluon model to compute the 
mixing of the scalar quarkonium and glueball states related to the observed /o resonances 
above 1 GeV. The model is contained in the interaction Lagrangian 

^ — -^g ~l~ ~l~ ^mix i (l) 

where Cq describes the quark sector, Cg the scalar glueball and Cmix introduces the mixing 
between these configurations. In the following we present details for the effective interaction 
Lagrangian describing the coupling between scalar mesons and their constituents. 



A. Quarkonia sector 

The quark-quark interaction for the isoscalar nn = + dd) and ss meson states 

above 1 GeV is set up by the isospin symmetric Lagrangian 

A = §^ (J'ni^) + (2) 

where the scalar quark currents are given by the non local expressions: 

Ux) = ^Jd'y {n{x + y/2)n{x - y/2)) ^{y) with n{x) = ( j (3) 



and 



Js{x) = Jd^y{s{x + y/2)s{x-y/2))^y). (4) 



The Lagrangian (2) is quartic in the quark field as in the NJL model [19-21], but with 
a non-local separable interaction as in [24]. The function represents the non-local 

interaction vertex and is related to the scalar part of the Bethe-Salpeter amplitude. In 
momentum space we have the Fourier transform of the vertex function, which 

characterizes the finite size of hadrons, is parameterized by a Gaussian with 

Hql)^exp[-ql/A% (5) 

where qe is the Euclidean momentum. This particular choice for the vertex function pre- 
serves covariance and wasjised in previous studies of light and heavy hadron system [22-25]. 
Any covariant choice for $ is appropriate as long as it falls off sufficiently fast to render the 
resulting Feynman diagrams ultraviolet finite. The size parameter A will be varied within 
reasonable range, checking the dependence of the results on it. 

The Lagrangian (2) represents the isoscalar sector of the more general quartic nonlocal 
S'[/y(3) flavour-symmetric Lagrangian [24] 
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c^'^j:{rf (6) 



0=0 



where J"(x) = / d'^yq{x + y/2)\^q{x — y/2)^{y'^) with g* = -u, d, s ^ , A" are the Gell- 
Mann matrices with A" = J2/31. In the isoscalar sector, i.e. considering the a = and 



a = 8 components, we are left with (2), which we expressed in terms of the yl/2{uu + dd) 
and ss configurations. In this work we do not consider the axial 5*^74(3) transformations, 
i.e. we do not relate the scalar quarkonia nonet to the pseudoscalar one. Although this 
would be a desirable feature, the connection of these nonets is not easy to establish. Such 
an approach with a chirally invariant interaction is for example pursued in the NJL model, 
but the predicted scalar masses are below 1.5 GeV [19,20] in conflict with the observed 
experimental resonances. In the following we restrict the study to the scalar sector, more 
precisely on the scalar-isoscalar states, which in turn can mix with the glueball configuration. 

The coupling strength of the current-current interaction is denoted by Kjsf and is in turn 
related to the meson masses. The scalar meson masses Mn and Ms are deduced from the 
poles of the T-matrix and are given by the solutions of the Bethe-Salpeter equations: 

The mass operator E;v(5)(p^), where p is the meson momentum, is deduced from the quark 
loop diagram of Fig. 1 and given by 

Siv(5)(/) = -tNcJ [Snis){q + p/2)S^is){q-pm] ^'{q'); (8) 

where A^c = 3 is the number of colors and Sn{s) is the quark propagator. A Wick rotation 
is then apphed in order to calculate Ejv(5)(p^). As evident from (7), Mjv and Ms are not 

independent. Once the size parameter A is chosen and using the bare nonstrange meson 
mass Mjv as an additional input, the coupling constant Kj^ and the bare strange meson 
mass Ms are fixed. 

We now turn to the discussion of the quark propagator. The general form [26] is 

Considering the low energy limit we have Z{p^) ~constant and m{j?) ~ m*, where m* is 
the effective or constituent quark mass. These masses are typically in the range of 0.25 to 
0.45 GeV for the u{d) quarks and 0.5 to 0.7 GeV for the s fiavour. The low energy limit 
for the quark propagator is too naive for our purposes, since it leads to poles in the mass 
operator for meson masses of about 1.5 GeV . In the following we consider two possible 
ways to avoid these infinities. The mass function for n = u,d fiavor can be decomposed as 
fT^nip"^) = '"^in+^n(p^) and the contribution of o"„ is assumed to be replaced by a large average 
value. To avoid poles, that is the unphysical decay of the quark-antiquark meson into two 
quarks, the averaged mass function is chosen as (mn(p^)) = m* + (fTn(p^)) = //n > 0.86 GeV 
(here for the u- and d -flavour), where 2 • is a lower limit set by the mass of the /o(1710) 
resonance. We then have 
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We also introduce an analogous parameter for the s flavour, which for jig ^ fin contains 
flavor symmetry violation. Above parametrization of the quark propagator is the simplest 
choice; although it neglects the dependence of the quark mass, the approximation by a 
free quark propagator with a large effective mass allows to test the approach and leads to 
considerable simplifications concerning the technical evaluation. A similar quark propagator 
was also used in Ref. [16], where the decay of scalar mesons into two photons was analyzed. 
We also consider a quark propagator which is described by an entire function: 

Siip) = (l - exp {(3 - n^))) . (11) 

This parametrization has been used both in the study of meson and baryon properties [27,28] 
and serves as one possible way to model confinement; the factor multiplying the free quark 
propagator removes the pole and on-shell qq creation is avoided. The effective quark masses 
rrii (with i^n,s) are taken from [29] with m„ = 0.462 GeV, rus = 0.657 GeV. In Ref. [29] the 
effective quark masses are calculated within a Dyson-Schwinger approach and, as functions 
of the Euclidean momentum, display an almost constant behavior up to high values. The 
parameter f3 is constrained from below to generate a behavior like i/ — rrii) for small (and 
for euclidean) momenta; (3 is also constrained from above, requiring that the propagator 
does not diverge for momenta up to ~ 2 GeV. 

The preceding discussion concerning the set up of the quark interaction Lagrangian 
(2) and the resulting generation of quark-antiquark meson states can also be alternatively 
described. By introducing the auxiliary meson fields A^(,7:) and S{x) {N for the ^J^(uu + dd) 
and S for the ss meson) previous procedure can be summarized by the Lagrangian [20]: 

4 = gNN{x)Jn{x)+gsS{x)Js{x) (12) 

with the condition that 

9n(S)-[ w-^ ■ (13) 

The last relation is known as the compositeness condition, originally discussed in [30] and ex- 
tensively used in the study of hadron properties [31]). The compositeness condition requires 
that the renormalization constant of the meson fields N{x) and S{x) is set to zero, hence 
physical meson states are exclusively described by the dressing with the constituent degrees 
of freedom. £^ contains the meson quark-antiquark vertex with a momentum-dependent, 
non-local vertex function. The resulting mass operator is completely analogous to the bub- 
ble diagram of Fig. 1. The Lagrangian (12) is very useful to calculate the decays of the 
bound state N{S); the coupUng constant gN{s) directly enters in the decay rate. 

B. Glueball sector 

We now extend the model to include the scalar glueball degree of freedom, described 
as a bound state of two constituent gluons. First we consider the scalar gluonic current 
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J = F;^(a;)F"'^'^(y) where F;^(a;) = d^Al{x) - d^Al^{x) + gfabcA^Al and A°(a;) is the gluon 
field. To set up a scalar glucball state we restrict the current to its minimal configuration, 
such that the glueball is described by a bound state of two constituent gluons. With this 
truncation we have J — F^^{x)F°-'^''{y) where F^y{x) = d^j^A'^^x) — duA'^^{x) is the abehan 
part of the gluonic field tensor. The effective, non-local interaction Lagrangian for the 
constituent gluons is then written as: 

A = ^Jli^l (14) 

where 

J,{x) = j d'y {f;^{x + yl2)F'^'^^{x - y/2)) (15) 

The non-perturbative features of the gluon dynamics are assumed to be taken into account by 
the coupling constant Kq. It should be stressed that the gluons introduced in this section 
are not the "background" gluons responsible for confinement [32,33], but two effective- 
constituent gluons forming the glueball. For simplicity we use the same vertex function as 
in the previous case, setting the glueball size equal to the one of the quarkonia states. The 
coupling constant Kq again is linked to the bare glueball mass by the pole-equation: 



1 

Eg(M2, 



where the mass operator Eg'(p^), as indicated in Fig. 2, is given as: 

Eg(/) = i2{Nl - 1) / [8(gi ■ q^)' + D{ql)D{ql)¥{q\ (17) 



where gi = q-\- p/2 and q2 — —q-\-p/2. D{q^) is the scalar part of the gluon propagator, 
which in the Landau gauge is [34] 

<. = ^"'(^M^-^)^(a- (18) 

As before, the formalism can alternatively be defined by a Lagrangian containing the 
scalar glueball field G{x) with: 

C'^ = gGG{x)J,{x), (19) 
where the coupling go is deduced from the compositeness condition 

The last two defining equations are equal to the ones used in [35]. 
For the gluon propagator we choose the free one [34,37] 

D{q") = (21) 
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where the effective mass rUg should be large, with 2 • nig larger than the bare glueball 
mass deduced from lattice simulations between 1.4-1.8 GeV. [6,38-40]. In the Wilson loop 
approach of Refs [33,36] a similar scalar glueball mass, around 1.58 GeV, is found. 

The generation of a constituent gluon mass is analogous to the effective quark mass, and 
one can relate its value to the gluon condensate [41]. Typical values for the effective gluon 
mass are in the range of 0.6-1.2 GeV. Lattice simulations give a value of 0.6-0.7 GeV [42]; 
in the QCD formulation of [43] a higher value of about 1.2 GeV is deduced. A similar value 
is also found in [44] for the off-shell gluon mass in the maximal abelian gauge. The gauge 
(in) dependence of the effective gluon mass is still an open issue [37,45]. For our purposes 
it is sufficient to take a constituent gluon mass within the range of 0.6-1.2 GeV; we choose 
an intermediate value of 0.9 GeV, which is in accord with the effective constituent gluon 
mass found in the study of gluon-dynamics in [46] . In the quarkonia sector we use a quartic 
interaction Lagrangian; this means that the background gluons, responsible for the string 
tension, i.e. for the attraction among the q — q pair, are described by a constant propagator. 
The gluons, exchanged among the quark- antiquark pair, are "soft"; the gluon propagator of 
eq (21) becomes a constant for small momenta, thus being in accord with this assumption. 
The exchange of non-perturbative gluons is also responsible for the appearance of a vertex 
function, i.e. for the finite size of the quarkonia mesons. In the current model calculation 
the vertex function is simply parametrized by (5). 

These arguments are also vahd for the soft background gluons exchanged by the two 
constituent gluons forming the glueball. On the other hand, when evaluating the glueball 
mass-operator (Eq. (17) and Fig. 2), the two constituent gluons can also have large timelike 
momenta because of the large glueball mass (~ 1.5 GeV). In this case the use of a constant 
for the gluon propagator is not adequate. This is why we use the form of Eq (21) with an 
effective mass. 



We also have to introduce a term which generates the mixing between the bare glueball 
and the quarkonia states. Using the flavour blindness hypothesis, which states that the 
scalar gluonic current only couples to the flavor singlet quark-antiquark combination, we 
have: 



where no direct coupling between the scalar quark currents J„ and Jg has been taken into 
account. In Eq. (22) the two constituent gluons forming the glueball can transform into a 
constituent quark-antiquark couple forming the quarkonia meson. It should be stressed that 
Eq. (22) does not describe the fundamental quark-gluon interaction, but an effective coupling 
of the scalar glueball current with the scalar quarkonia ones, suitable for the description of 
the mixing among these configurations. 

Following the arguments of [1] a direct mixing of the bare quarkonia states is a higher 
order perturbation in the strong coupling eigenstates and is neglected in the following. 

Due to the mixing term physical states are linear combinations of the bare quark- 
antiquark and glueball states. Implications of this scenario will be studied in the next 
section. 



C. Mixing sector 




(22) 
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III. MODEL IMPLICATIONS 



A. Masses of the mixed states 



We now consider the total Lagrangian containing the interaction terms of (2), (16) and 
(22) with 



(23) 



Due to the mixing term Cmix the Bethe-Salpeter equations (7) and (16) are not valid any- 
more. In fact, the T-matrix is modified by Cmix and takes the form: 



T 



-{l-K- T?)'^K = -{K-^ - S)-\ 



(24) 



with the non-diagonal coupling matrix 



K 



/ V2Ksg 

V^KsG Kg 
\0 KsG 



\ 

KsG 



(25) 



and the diagonal mass operator 




(26) 



The masses of the mixed states are obtained from the zeros of the determinant in the 
denominator of the T matrix with Det [1 — Kiy\= 0. In the limiting case of Ksg = the 

determinant reduces to (1 — Kn^m{p^)){1 — Ks^s{p^)){.^ — Kg^g{p^)) = 0, which results 
in the defining equations of the bare masses without mixing (cqs. (7) and (16)). 

For Ksg 7^ we have mixing, where, starting from the unrotated masses (Mat, Mg, Ms) , 
we end up with the mixed states (A^', G', S') of masses (M^v', Mg', Ms')-, which we interpret 
as the physical resonances /o(1370), /o(1500) and /o(1710). Quantitative predictions in the 
three-state mixing schemes strongly depend on the assumed level ordering of the bare states 
before mixing. In previous works [1,8] essentially two schemes were considered: < 
Mg < Ms and M^ < Ms < Mg- In phenomenological studies latter level ordering seems to 
be excluded, when analyzing the hadronic two-body decay modes of the /o states [7,9]. In 
the current work we will refer to the first of these two possibilities, where the bare glueball 
mass is centered between the quarkonia states before mixing. We also compare our pole 
equation Det [1 — KTj\= to other approaches. But first we discuss further consequences of 

the covariant non-local approach, such as the meson-constituent coupling constants, which 
are crucial in the calculation of the decays of the mixed states. 
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B. Resulting meson-constituent coupling constants 



As a result of the mixing the rotated fields A^', G' , S' couple to the nn, gg and ss 
configurations. The leading contribution to the T-matrix (24) in the limit ~ Mf with 
i ^ N',G',S' is given by the pole at = Mf with 

rpab 9i 9i f07\ 

where a,b = nn, gg and ss refer to the constituent components. A similar expression for the 
T matrix is also given in [19,20], where the rj — r]' mixing was analyzed in the context of the 
NJL model. The constant g^^^, for example, represents the coupling of the mixed state G' 
to the quark- antiquark configuration of fiavor u and d. The modulus of the nine coupling 
constants is then given by 

which we solve numerically. In Appendix A, where we consider the mixing of two fields , we 
also give an explicit expression for this quantity. The expression (28) is the generalization 
of the compositeness condition (eqs. (13) and (20)). 

In the limit Ksg the T-matrix is diagonal and the coupling constants gQ,, for 
instance, become 

0, (29) 



nn 

9g' 


— 9g' 


\9g'\ 


— 9g 



Qp2 



(30) 



p^=Ml,=M% 



where the last equation is the glueball compositeness condition (20). In this case there is no 
mixing and the glueball couples only to gluons. 

We still have to discuss the sign of the coupling constants. By convention we chose (7™, 
gQt and g'gt as positive numbers. The sign is determined from the off-diagonal elements of 
the T-matrix: 

gt = szgn{<4) \J (31) 

where 

Having obtained these coupling constants we again can write an effective interaction 
Lagrangian for the G' field 

^'g' — 9g'^' Jn + 9g'^' '^9 + 9g'^' Jsi (33) 
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where the couphng between the mixed states and the constituent configurations are made 
exphcit. This Lagrangian allows to calculate the decay of G' as we will see explicitly for the 
two-photon decay; the coupling constants directly enter in the decay rate of the state. One 
then has completely analogous expressions for the fields A^' and S': 

C!^, = ^piV' J„ + g'^^N'Jy + (f^,N'J,, (34) 
Cls' = g'pS'Jn + g'SS'J, + ^iS'Js. (35) 

The strength of the coupling is directly connected to the mixing strength, the explicit 
mixing matrix is discussed in the following. 



C. Mixing matrix M 



When dealing with elementary scalar particles, and not with composite ones, mixing 
between the bare meson fields [N ^ G, S) can be expressed by the Lagrangian 



+fGS + V2frGN, 



(36) 



where / and r are mixing parameters (r 7^ 1 takes into account breaking of the fiavour 
blindness hypothesis). As before, no direct mixing between A^" and S has been taken into 
account. In this simple case one has to diagonalize 



K-G 



( -Ml V2fr \ 

V2fr -Ml f 
VO / -MlJ 



(37) 



where the mixing matrix Mk-g C 50(3) is obtained from the condition 
Mk-g ■ Tk-g ■ M],_G = Dk-g = Diag[-Ml,, -Ml,, -M|,], 



(38) 



where the primed Klein-Gordon masses refer to the mixed states. The matrix Mk-g — 
Mif_Q C SO{3) (i = A^', G', S' and a ^ N,G, S) connects the unmixed and the mixed 
states by 



(39) 



In the following we want to determine an analogous matrix in our Bethe-Salpeter ap- 
proach. We consider \G') which can be written as a superposition of the bare states: 







(\N)\ 


\G') 


^Mk-g 


\G) 


V 1^') ) 




\\S) 1 



G') = M^''^ lAT) + M^''^ \G) + M^'-^ \S) 



g'g 



tG',S 



(40) 



where M^''^ , for example, is the admixture of |A^) = uu + dd^ to the mixed state \G'). 
Unlike in the Klein-Gordon case care should be taken with such an expression, since here the 
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bare states |A^) , IG) and \S) are not well defined. They correspond to the Bethe-Salpeter 
solutions in the case of zero mixing, but, when including mixing, they are not normalized 
vectors of the Hilbert space anymore. To obtain a corresponding expression as in the Klein- 
Gordon case, we proceed as follows: the coupling constant Qq} (evaluated in the previous 
section) is related to qn (the couphng constant of the bare state N to nn in the case of no 
mixing) by the relation 

^g'/ = Af«'^^?^(p2 = M^,), (41) 

where QNip^ = Mq,) is determined from (13) evaluated at the physical mass with p'^ = Mq,. 
Eq. (41) states that the nn to G' coupling equals the nn admixture in G' (which is the 
matrix element we want to determine) times the N to nn coupling evaluated on 

mass-shell of G' . 

Using (13) and (31) we obtain for M^''^ : 

-1/2 

(42) 

Generalizing this result for a generic component of the mixing matrix M we have 

-1/2 

(43) 

with i — N', G', S' and a — N = nn, G = gg, S = ss. Last expression, while not being 
based on a rigorous derivation, can be regarded as a definition of the mixing matrix in a 
Bethe-Salpeter approach. We explicitly show in Appendix A for the reduced problem of two 
mixed fields that this definition of the mixing matrix is the analogous one the Klein-Gordon 
case. 

Furthermore, the components of the mixed state |i) are correctly normalized with 

1 = (M'^f + {M'^f + (M'^)2 (44) 

as verified both numerically and analytically (Appendix A). This is a further confirmation 
of the consistency of the definition (43). While the rows of M are properly normalized, 
this does not hold for the respective columns of M because of the p"^ dependence of the T 
matrix and of the mass operators. The rows are evaluated at different on-shell values of 
values with = Mf. This imphes that the matrix M is not orthogonal, but here we 
will demonstrate that deviations from orthogonality are small for both choices of the quark 
propagator. 

D. Comparison with other mixing schemes 

We again consider the Klein-Gordon case and its mass equation Det[p'^ + Tk-g] = 
which reads: 

(/ - O ip' - Ml) ip' - Ml) - f ((p2 - M^) + 2r' ip' -Ml))^ 0. (45) 
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9g' 



gr,{p' = Ml,) 



= sign{aQ}) 




M' 



ga{p' = Mi) 



sign{a1) 



Qp2 



Qp2 



In our approach the mass equation Det[T ^] — (see Eq. (24)) is written out as: 

(1 - ir^S^(p2))(l - KcJ^Gip'm - KsMp'))- 



+2Eg(p')E5(p')(1 - KsMp'))} = (46) 

In order to compare the last expression with (45) we introduce the functions rjaip^) as 

(1 - K^Eaip')) = ip' - M'Mp') (47) 

for each a = N,G, S. The rjaip^) do not contain poles for = since (1 — KaEaip^ — 
Ml)) = (see (16) and (7)). When substituting (47) in (46) we get: 

{p'-Ml){p'-Ml){p'-Ml)- 

2 SG(p^)Sg(p2) / 2 ..2^ I o Sjv(p') r]s{p^) , o ,.2^^ r. /.oN 

vg{p^)vs{pV \ m{p ) ^s{pV J 

Comparing (48) with (45) we deduce that the mixing parameters / and r of the Klein- 
Gordon approach become p^-dependent functions in the case of composite scalar fields. In 
particular we make the following identification 

p2, .^ _ r.2 ^g{p')^s{p'' 



m{p'^)vs{p^) ' 



riN{p')^s{p^r ^ ^ 

It should be noted that the composite approach generates a value of r 7^ 1 which reflects 
the deviation from the flavor blindness hypothesis. Although we set up the quark-gluon 
interaction assuming this hypothesis, on the composite hadronic level we obtain a breaking 
of this symmetry due to the flavor dependence of the quark propagators. As we will see in 
our numerical evaluation (sec section IV and Figs. 5 and 6) the functions /(p^) and r(p^) 
vary slowly in the momentum range of interest, thus justifying a posteriori the Klein- Gordon 
approach. 

Replacing in (37) the constants / and r by the running functions f{p^) and r(p^) we 
have 



/ -Ml V2f{pMp') \ 

- TK-aip") = ^/^f{pW) -Ml f{p^) 

\ m -Ml ) 



(50) 



where the Klein-Gordon mass equation Det\p'^ + Tk-g{p'^)] = coincides with the pole equa- 
tion Det[T~^] = by construction. When diagonahzing Tk-g{p^) we obtain a momentum- 
dependent transition matrix Mk-g{p^) C (SO (3) (for each p^), which cannot be directly 
interpreted as the mixing matrix. To extract the composition of a mixed state one has to 



12 



consider its on-shell mass value. For example, the mixing coefficients for the diagonal state 
\G') arc obtained from the second row of Mx-cip^ — ^g') ^^c corresponding on-shell 
value. We then can define a mixing matrix M' (the prime serves to distinguish it from M 



defined in the previous subsection) where we have M'^'^' 
and in general we get: 



M 



G',a 
K-G 



(/ = M2,) {a = N,G,S) 



(51) 



Above procedure arises from the analogy to the Klein-Gordon case and indicates that M' 
is a natural generalization of Mk-g, provided that we evaluate the rotated states at their 
corresponding on-shell mass value. Again, due to the p^-dependence and the on-shell eval- 
uation the mixing matrix M' is not orthogonal, where in the limit of /(p^) = const and 
r(p^) = const one has M' — Mk-g as desired. The numerical results of M (43) and M' 
(51) are very similar, as we will show in the next section. The analytical argument leading 
to M ~ M' are presented in Appendix A. Because of the weak dependence of the mixing 
functions /(p^) and r{p'^), the mixing matrix M (and M' ~ M) will be almost orthogonal, 
as shown in section IV. 

In the three-state mixing scheme many phenomenological approaches [1,7,9] refer to an 
interaction matrix, which is diagonalized linearly with respect to the masses: 



Tc 



QM 



f -Mn V2f*r* \ 

x/2rr* -mg r 
\o r -Ms) 



(52) 



Here we introduce the parameters /* and r* to distinguish them from those in (45). The 
resulting equation for the masses of the mixed states is then: 



Mjv)(i 



MG){Jp^-Ms)-f 



*2 



Mn) + 2r 



*2/ 



Ms) = 0. (53) 



To relate our approach to the linear mass case with parameters /* and r* we start from (48) 
and decompose (p^ — Mf) — ( VP^ — M^) ( VP^ + Mi). Then we have 




p2 - Mjv)(Vp' - Mg)(Vp' - Ms) 
p^ — Mn) + 2r (p )- ' 



fHir) 



iVp" + Ms){Vp' + Mg) 
f-Ms)] =0. 



[Vp' + Mg) 

Comparing (54) to (53) we find the "running" behavior for /* and r* with 



(54) 



y + Ms)iV^ + MG) 



„*2 



(V^ + Mg) 
(VP + Mg)' 



(55) 



where extra-factors occur when relating the quadratic to the linear mass case. The consis- 
tent hmit of our approach is actually the Klein-Gordon case, but Eq. (55) yields a direct 
comparison to the linear mass case, where many phenomenological approaches are working 



m. 
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E. Decay into two photons 



In this work we also analyze the two-photon decay rates of the scalar isoscalar /q mesons 
considered. The two-photon decay width constitutes a crucial test to analyze the charge 
content of the scalar mesons [9,47,48] . The glueball component does not couple directly to 
the two-photon state and leads to a suppression of the decay width when present in the 
mixed state. 

For the decay of the N' = /o(1370) meson by its uu constituents into two photons we 
have to consider the triangle diagram of Fig. 3. Analogous diagrams occur for the d and s 
flavors in the loop. 

When evaluating the two-photon decay we will consider the first choice, that is the free 
form of the quark propagators, which does not imply a modification of the QED Ward 
identity, thus considerably simplifying the problem. Since we deal with a non-local theory 
as introduced by the vertex functions, care has to be taken concerning local gauge invariance. 
The triangle diagram amphtude of Fig. 3 

MZangle = / " Tr[SMr Sn{p2)r Sn{p3mq% (56) 

where ki and k2 are the photon momenta and pi — k + ki, p2 — k, — k — k2, is not gauge 
invariant. 

The trace reads 

= 4/x„ {{k-(k^ - {ki ■ k2)gn) + ^l^ni^k^k^ - g^'ik' - i^D) (57) 

where we have omitted the terms proportional to k^ and k2 because they do not contribute 
to the decay. The first part of the trace is clearly gauge invariant, while the second one is 
not. Gauging in addition the non-local interaction term (for the procedure see [22,23,49]) 
gives rise to the two extra diagrams of Fig. 4, where the photons can now couple directly 
to the non-local vertex. Each separate diagram is not gauge invariant, but the physical 
amphtude, which is the sum of the three diagrams, is gauge invariant [22] . We can find the 
gauge invariant contribution of the triangle diagram by considering the shift 

Y - 7^^,., =Y- (58) 
7^^ - i1m = 7^^ - (59) 

K2 

The term 

^r.U^e,± = • / (0 ■ ^^['^"(^^i)^'>^i'^"(^^2)^^.^.'^-(^^3)]^(g') (60) 

is then by construction gauge invariant. As shown in [22], which we refer to for a careful 
analysis of this issue, this can be done for each diagram and the total amplitude is 

^total — ^triangle + ^bubble + -^tadpole — ^triangle,± + ^bubble,± + ^tadpole,± (^1) 
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where the gauge breaking terms cancel. 

As already discussed in [22] and also checked for our case numerically, the gauge invariant 
part of the triangle diagram -/W^riansZe.x represents the by far dominant contribution to the 
two-photon decay width. For the case of a Gaussian vertex function the bubble and tadpole 
diagrams of Fig. 4 are further suppressed in the order of 10~^. 

We finally give the analytic expression for M^'^^^^^^g 

Mtr-angle,± = {^^'2 " (^1 " ^2)9^ ' I, (62) 

where I — Ii + I2 is calculated from (60) with 

The term Ji also occurs in the neutral pion decay into two photons, whereas the additional 
term I2 turns out to have opposite sign to Ji and tends to lower the decay rate. 

With above results we finally obtain for the two-photon decay width of the mixed scalar 
meson state hke N': 



(65) 



where Nc = 3 is the number of colors and I^^^ is obtained from I by replacing fin with Hs 
in the quark propagator. Only the coupling constants g^^f and g'f^, of the mixed field A^' 
contribute to the decay rate, while g'j^, does not because gluons do not couple directly to 
photons. The extra-factor 2 in front of the coupling constants comes from the exchange 
diagram. Analogous expressions follow for the decay rates of the other two resonances 5" 
and G'. 



IV. NUMERICAL RESULTS AND DISCUSSION 



First we proceed to constrain the parameters of the model. The mass scale of the non- 
strange quarkonia in the scalar meson nonet is set by the physical state ao(1450). We 
therefore assume that the mass of the unmixed N meson state is close to this value. Using 
Mtv = 1.377 GeV as an input, which is the value deduced in the phenomenological analysis 
of Ref. [9], the mass of the bare S meson state and the elementary quark coupling constant 
A'at are fixed from Eq.(7) once the cut-off A is specified. The bare scalar glueball mass is 
predicted from lattice calculations in the range of 1611 ± 30 ± 160 MeV [6]. The analysis 
of [9] prefers a value of 1.443 GeV, which we use as an additional input. The bare glueball 
mass fixes in turn the gluonic coupling Kq of Eq. (16). 

We now consider separately the phenomenological consequences of the two proposed 
choices for the quark propagators. 
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1. Case 1, free quark propagator 



a. Mass spectrum: As discussed in Section II we consider large effective quark masses 
in order to avoid poles in the integration. For the quark mass we choose the threshold value 

— 0.86 GeV to prevent unphysical on-shell production of a quark- antiquark pair. 

For the vertex function wc choose a cut-off of A = 1.5 GeV comparable to the values 
chosen in the analysis of [22] for the hght meson sector. We will subsequently vary A 
within reasonable range, in order to check the dependence of the results on the specific 
value. Actually, we find a remarkable stability of the results under changes of A as discussed 
further on. The mixing coefficient Kga and the effective strange quark mass iXg in the quark 
propagator are fixed by the experimental masses Mq' = 1507 ± 5 MeV and Aigi = 1713 ± 7 
MeV. We then can compare Ksa with phenomenological approaches and with estimates 
from the lattice following the expressions of (49) and (55). Here we do not use the mass of the 
resonance N' = /o(1370) as a further constraint, since it is rather broad and ill-determined. 
For the optimal choice of /ig — 0.989 GeV and Ksg — 0.55 GeV~^ we obtain for the masses 
of the physical N' and the bare S states: 

Ms = 1.696 GeV, Mn' = 1.297 GeV (66) 

The obtained quark "mass difference" — t^n = 129 MeV is close to the upper limit of the 
current quark mass difference of about 130 McV. With these quark mass values the bare 
mass level scheme Mn < Mq < Ms comes out naturally. A reversal of the bare scheme 
with Miv < Ms < Mq would require a rather small mass difference — IJ'n in conflict with 
phenomenology. 

b. Mixing matrix: The mixing matrix M linking the rotated and the unrotated states 

is: 

0.80 0.59 0.10 \ (\N)\ 

-0.60 0.76 0.26 1^) (67) 
0.08 -0.27 0.96 / V i-^) / 

showing a similar pattern as in the phenomenological work of [9]: the center state \G'), 
identified with the /o(1500), has a dominant glueball component and the quark components 
\N) and IS*) which are out of phase. Latter effect causes destructive interference for the KK 
decay mode consistent with experimental observation [1]. On the other hand |A^') and |iS") 
show a tendency to be dominated by the nn and ss constituent components as also deduced 
in [7,9]. 

As already indicated, the mixing matrix of (67) is not orthogonal. The deviation from 
orthogonality is displayed by M ■ (which is just the identity in the Klein-Gordon limit): 




M-M* 



/I -0.0028 -0.0019 \ 

-0.0028 1 -0.0074 
V -0.0019 -0.0074 1 / 



(68) 



where the off-diagonal elements turn out to be very small. This in turn implies that the 
Klein-Gordon limit is fully appropriate to set up the mixing scheme of the scalar meson 
states. Furthermore, the unities on the diagonal are in accord with (44). 
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In Section III.C we have introduced the mixing matrix M'. The numerical evaluation 
yields 



M' = 



/0.79 0.60 0.09 \ 

-0.61 0.75 0.24 
V 0.08 -0.26 0.96 / 



(69) 



where the result is nearly identical to the one for M. The reasons for this coincidence are 
given in Appendix A. 

The pattern observed in the mixing matrix only weakly depends on the particular choice 
of A. Similarly, a change in the value of /x„ up to 1.1 GeV alters the mixing pattern only 
within 5 %. The weak dependence of the results on A and /x„ is explicitly indicated in 
Appendix B. 

c. Coupling constants: In the quantum mechanical approach the amplitude for the 
decay of the mixed state \i) (with i = N' , G' , S') into two pions is fed by the > component 
in the strong coupling limit [1] . Hence the 27r decay amplitude is proportional to Af defined 
in Eq. (43). In this limit wc obtain for the ratio of 2n decay amplitudes \An/^2it/Ag'^27t\ 
a |m^''^/M^''^| = 0.80/0.60 = 1.33 (in the preferred solution of [9] one has for this 
ratio 1.13). In our full approach the decay amplitude is related to the respective coupling 



„nn /„' 

9n'I 9, 



1.64, which shows a 



constant, hence we obtain for the ratio |AAr'-»27r/^G'-^27r| o; 
clear enhancement. This might explain why the state N' = /o(1370) is considerably broader 
than the G' = /o(1500). The strong deviations of mixing coefficients and coupling constants 
is a non- negligible effect, which can be traced to the fact that the coupling constants Qaiv^) = 
(dHaip^) / dp^) with a = N,G, S are momentum dependent. If Qaip^) = const, we would 
recover the quantum mechanical limit of /M^ — g'^/g(p (and similarly for all the 

other ratios). Again, the relation g^ / g^l' > M^''^/M'^''^ is stable when changing the 
parameters A and fin- 

The complete set of resulting coupling constants is summarized as 



„nn . „nn . „nn 
9n' ■ 9g' ■ 9s' 

9n' ■ 9g> ■ 9f' = 

9n' ■ 9g' ■ 9s'' -■ 



-- 7.35 : -4.48 : 0.40 \ 
1.20 : 2.69 : 8.10 
: 0.91 : 0.95 : -0.23 y 



(70) 



These results indicated here are of course model dependent, but they point to an inter- 
esting aspect in the comparison with other studies: the ratio of coupling constants, relevant 
for the strong two-body decay modes, can vary rather sensibly from the ratios of mixing 
matrix elements. This deviation is a consequence of the bound-state nature of the scalar 
mesons in a covariant framework. 

It would of course be interesting to calculate the decays into two pseudoscalar mesons 
by loop diagrams directly. But this would require a consistent knowledge of the quark 
propagators and a careful study of the pseudoscalar meson-quark vertex functions to treat 
both the scalar mesons above 1 GeV and the light pseudoscalar mesons in a unified way. 

d. Running mixing functions: As already explained in the section devoted to the com- 
parison with other mixing schemes, here we obtain a momentum dependent mixing strength. 
We developed two formulations, /(p^) of Eq. (49) and f*{p^) of Eq. (55), in order to com- 
pare to the Klein- Gordon case and the quantum mechanical linear mass limit. Our results 
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for p^-dependence of the mixing strength is summarized in Fig. 5. The running function 
f*(jp') varies in the range between 60 and 64 AdeV for the values of interest. Our result 
should be compared to the value of 85 ± 10 MeV obtained in [9]. Lattice calculations in 
the quenched approximation obtain /* = 43 ± SlMeV [8,38] with a large uncertainty but 
of similar magnitude. Other quantum mechanical studies use fitted mixing parameters of 
/* = 77MeV [8], /* = 64 ± 13Mey [38] and /* = 80 MeV [7]. Since f*{p'^) does not vary 
drastically in the region considered, this weak dependence explains why the mixing matrix 
M (and analogously M') is "almost" orthogonal. In the effective QCD approach of Ref. 
[32] a smaller quarkonia-glueball mixing is found; however, strong mixing is obtained when 
including an intermediate hybrid state appearing in the same mass region as the glueball 
and quarkonia ones. 

Another characteristics of the non-local covariant approach is the dynamical generation 
of flavour blindness breaking with r > 1. The values of the matched "running" functions 
r(p^) and r*(p^) (Fig. 6) are by 10-20 percent larger than unity. The lattice result of 
r* — 1.20 ± 0.07 [8] is in rather good agreement with our evaluation. 

The characteristics of the mixing parameters are again rather stable when changing the 
parameters. 

e. Two-photon decay widths: We first consider the decay widths of the bare states 
and S when mixing is neglected. In this case the coupling constants N — nn and S — 'ss are 
given by (13). We obtain 



0.821 keV, Ts^2'y = 0.08 keV. (71) 



As explained previously, the amplitude of the triangle diagram which gives the dominant 
contribution to the decay is proportional to (Ji + 12). The extra term I2, absent in the 
neutral pion case, has the opposite sign to Ii and the ratio Ih/h] grows with increasing 
mass of the resonance. In the decay of the bare state N the term I2 lowers the decay rate 
by a factor of 5.22. 

For the two-photon decays of the physical states, where the coupling constants of the 
mixed states are used, we get: 

r;v'-.27 = 0.453 keV, rG'-.27 = 0.273 KeV, rs'^2'y = 0.125 keV, (72) 

resulting in the ratios of 

rjv'^27/r5'^27 = 3.62, rG"^27/r5'^27 = 2.19. (73) 

which are similar to the ones obtained in [9] with rAr/_,27/r5/^27 = 3.56 and rG"-»27/rs''^27 = 
2.36. Current experimental upper limits for G' = /o(1500) and for S' = /o(1710) are [50]: 



/o (1500)^27- 



/o(1500)-^7r7r 
r/o(1500)tot 



< 0.46 keV, (74) 



r/o(mo)^27 < 0.110 keV . (75) 

i /o(1710)tot 

Multiplying our theoretical results by the experimental ratios [50] = 0.454 ± 

0.104 and ^/odno^KK ^ Q_s8t^-09 g^^. 
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0.124 ± 0.028 keV (76) 

0.04751^23? keV, (77) 

in accord with the experimental upper limits. 

The experimental two-photon decay width of the scalar resonance /o(1370) has been 
seen; originally on [51] two values were indicated, i.e. 3.8 ± 1.5 keV and 5.4 ± 2.3 keV. 
However, it is not clear if the two-photon signal comes from the /o(1370) or from the high 
mass end of the broad /o(400 — 1200). The PDG currently [50,52] seems to favor this last 
possibility, but it states in a footnote that this data could also be vahd for the /o(1370). We 
therefore interpret the two experimental values as an upper limit for the two-photon decay 
width of the /o(1370). The result for r7v'=/o (1370)^27 is an order of magnitude smaller than 
these upper limits. A precise experimental determination of the two-photon decay values of 
these scalar states would clearly help in understanding their structure. 

For what concerns the cut-off dependence of the decay rates, we note that an increase 
in the cut-off leads to a weak increase of the decay widths as indicated in Appendix B, but 
the ratios remain stable. 



(r/o(1500)^27)t/ieori/ 



r/o(1500)-»7r7r 
r/o(1500)tot 

foil710)^2-y )theory—^ 

i /o(1710)tot 



2. Case 2, entire function 

In the following we summarize our results for the quark propagator of Eq. (11), described 
by an entire function modelling confinement. 

a. Mass spectrum: We fix the parameters in the same fashion as for the previous case 
(that is to Mn = 1.377 GeV, Mq = 1.443 GeV, = 1.507 GeVand Ms' = 1.713 GeV), 
but this time we have the parameter /3 of the quark propagator of Eq. (11) together with 
the mixing strength Ksg- For an identical cut-off with A = 1.5 GeV we get: 

Ms = 1.698 GeV, Mn' = 1.297 GeV (78) 

obtained for the fit values of /3 = 2.44 GeV'^ and Ksg = 0.284 GeV'^. The two results 
of Eq. (78) are essentially identical to the previous case of Eq. (66). Also in this case the 

reversed bare level ordering with Mjv < Ms < Mq is disfavored, requiring a decrease of (3 
of the order of 10~^. This is in contrast to the requirement that the propagator behaves as 
a free in the limit of small Minkowsky or Euclidean momenta. 

h. Mixing matrix: The mixing matrix M linking the rotated and the unrotated states 
is practically unchanged when compared to case 1: 

N')\ / 0.80 0.59 0.10 \ /|iV)\ 

G') = -0.59 0.76 0.25 |G) , (79) 
S') j \ 0.07 -0.27 0.96 / V l-^) / 

where M ■ M* is very close to the identity matrix. 
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c. Coupling constants: The results concerning the admixture and couphng constant 
ratios is analogous to the previous case. For the ratio of mixing amplitudes we get 



N 



0.79/0.59 = 1.36, while for the coupling constants we have g^^/g^ 



1.65, again almost unchanged with respect to case 1. The complete table reads: 



Un' ■ yg' ■ ys' 

„ss _ 

9s' - 



9n' ■ 9g' 



\ 



„99 

9n' 



■ 9g' 



= 3.82 : -2.32 : 0.22 
0.63 : 1.37 : 4.11 



= 0.91 : 0.95 : -0.23 



(80) 



d. Running mixing functions Results for fip^), f*{p'^), r[p'^) and r*(p^) are summarized 
in Fig. 6, with the same quantitative behavior as in Fig. 5. 

We conclude this section by noting that the quark propagator, modelling confinement, 
gives rise to very similar results as the free propagator with a large effective quark mass. 
Again, changes in the cutoff A do not alter the quahtative features of the results. This result 
is rather encouraging, since the model predictions considered here, seemingly do not depend 
on the particular choice of the quark propagator. 



V. CONCLUSIONS 



In this work we utilized a covariant constituent approach to analyze glueball-quarkonia 
mixing in the scalar meson sector above 1 GeV. We used simple forms for the quark 
and gluon propagators in order to avoid unphysical threshold production of quarks and 
gluons. Although quark and gluon propagators are directly accessible in lattice simulations 
in the Euclidean region, an extrapolation to the Minkowsky region, also needed here, is not 
straightforward. We therefore considered relatively simple choices of the propagators, which 
allowed us to point out some features of the mixing of covariant bound states. We tried to 
work out similarities and differences with phenomenological approaches, in particular with 
respect to the analysis of [9]. Although in a covariant approach the mixing matrix is in 
general not orthogonal, in the present case only small deviations from orthogonality are 
obtained. This in turn leads to a mixing pattern rather similar to that of Ref. [9]. The 
mixing matrix M has been introduced by knowledge of the coupling constants of the mixed 
and unmixed states taking into account their p'^ dependence. The resulting matrix M is 
analogous to the Klein-Gordon case as shown both numerically and in part analytically. 

Many properties we analyzed, such as the appearance of the "running" mixing parameters 
/(p^) and r(p^), are rather independent on the choice of the particular quark propagator. The 
numerical results for /(p^) and r(p^) are in qualitative accord with the lattice evaluations. 
We generate a dynamical breaking of the glueball flavor blindness corresponding to r slightly 
bigger than unity, which is directly connected to the isospin symmetry violation at the level 
of the quark propagators. Again, all these considerations do not depend on the choice of 
the two proposed quark propagator forms and on the employed parameter sets. Another 
interesting result of our approach is that the bare level ordering Mn < Mq < Ms is naturally 
favored for both propagator choices. 

As a further application we also evaluated the two-photon decay rates in the context of 
the mixing model. The predicted results are in accord with the present experimental upper 
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limits. The respective ratios are also in agreement with the phenomenological estimate of 
[9]. 
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APPENDIX A: MIXING MATRICES M AND M' 



In Section III.C we introduced the mixing matrix M (see eq. (43)) and then in Section 
III.D alternatively M' (see eq. (51)). For both mixing matrices, as indicated in Section IV, 
we find similar results. To demonstrate the close connection between M and M' we consider 
the reduced problem for the mixing of two fields, G and 5, thus leaving out the field N from 
the discussion. In this case we have the rotated fields G' and S' only, and for G' we have 
\G') = M^'^ \G) + M^'^ \S) . and E are now 2x2 matrices with: 



K 



Kg Ksg 
KsG Ks 



Sg 
S5 



Then we write the full expression for the matrix T ^: 



from which we get 



Ks 



1 



Det[K] 
Ksg 
Det{K] 



Kr, 



+ Eg 



K. 



Det[K] 

Kg I V, 
Det[K] ^ 



Det[T-^] 



Det[K] 
Ksg 



Det[K\ 



(Al) 



(A2) 



(A3) 



Det[K] Det[K] 

The coupling constant Qq, (31) is then explicitly given as 



gg 
9g' 



G 



p2^M2, p2^M2, \ Det[T-^ ^ Det[K] 



+ E5) = 



fdDet[T-^] 



-1 



G 



Det[K] 




(A4) 
(A5) 



The expressions for the other coupling constants follow from (31). Similarly, the explicit 
expression for M*^ is 



G'G 



'1 



99 
9g' 



G 



Qp2 



G' 



( dDet[T-^] 



G 



Det[K] 



+ Ec 



G 



(A6) 



G' 
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One has then similar relations for the other elements. 

First we show that the components of the state \G') — M'^ ^ \G)-\-M^ ^ \S) are correctly 
normalized. Using (A6) and (A3) we have 



G'5\2 



dDcf[T 



Det[K] 



Det[K] 



■ (AT) 



G' 



We come now to the equivalence of M and M' . Let us consider the ratio ^ = M^'^ /M^'^ 
which can be calculated explicitly from the basic definition (43) and form (A3) : 



G 



-Kg + DetjK] ■ Eg 

KsG 



dT^c/dp^ 
\ dEs/dp^ 



(A8) 



(A9) 



where the last term has been obtained making use of the equation Det[T~^] — 0, which 
holds for = Mq, . 

When we consider the analogous ratio calculated from the elements of M' (43) (in this 
case we just have the running parameter /(p^) from equation (49) and not r(p^)) we find: 



a 



M'G'S 



KsG \ \T.g{p^)^s{P^))^: 



(AlO) 



We want to prove that ^ '^^ (i.e. M ~ M') in the weak mixing hmit {Ksa small) and in 
the f{p^) = const hmit, for which M' = Mk-g. 

A small mixing strength Ksg implies that we can neglect the term proportional to K^q 
in the determinant of equation (A9). Introducing the quantity rjsip^) of Eq. (47) we find: 



Kg{MI - Ml,)rjs{p') 



K. 



sg 



^[dj:s(p')/dp')^^. 



(All) 



For ~ Mq the following approximation is valid: 



In fact, when Ksg is small, Mq/ is close to Mq. This allows us to write 



(A12) 



.^Gip' 



(A13) 
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Let us now consider the second limit, for which /(p^) = const in the region of interest, 
i.e. between Mq, and M|, in the 2 field case, and between Mfj, and A'/|, in the 3 field case. 
The condition /(p^) = const is satisfied if ^a{p^)/Va{p'^) = Cq, where is a constant for 
a — G,S. In this case one has M' — Mk-Gi which is orthogonal. We are then considering 
the orthogonal limit for M' . The condition T,a{p^) / rjaip^) — Ca for the case a — S implies 
the following form for T^sip'^) '■ 

This is of course an approximate form for 'Es{p'^) valid in the limit of a constant cg. Note 
that the condition E5(p2 = M|) - l/Kg = from Eq. (2) is fulfilled. With this form for 
T,s{p'^) we have 



Vsip'^ 



Vsip"^ 



(A15) 



which is valid in the interval where (A 14) is valid. 

Plugging the approximations (A13) and (A15) in (A9) we have indeed 



^ ^ ^, ^ Ml - M l, 

KsG 




Ml - Ml, 



(A16) 



thus having similar results for M and M' = Mk-g for the used approximations. As shown 
in section IV, in the three field mixing case, similar matrices are found. This is then true 
for all the parameters studied in our work and for both propagator choices. 

Note that the two conditions discussed in this appendix are satisfied: Ksg is small 
(corresponding to a small difference Mg, — M| ) and the function f{p^) is almost constant 
in the momentum interval between M^, and M|, (see Fig. 5 and Fig. 6). 



APPENDIX B: RESULTS FOR PARAMETER VARIATION IN THE CASE OF 
THE EFFECTIVE FREE QUARK PROPAGATOR 

To explicitly indicate the dependence on the cut-off value of the vertex function here we 
summarize our results for A = 2 GeV with the fit values = 0.985 GeV, Ksg = 0.122 
GeV^^. The masses of physical N' = /o(1370) and bare S state are Mjv = 1.287 GeV and 
Ms = 1.696 GeV. 

Mixing matrix: 

/0.79 0.61 0.11 \ 
M = -0.61 0.74 0.27 . (Bl) 
V 0.08 -0.28 0.95 / 

Set of coupling constants: 
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/ 






■ 9s' 


= 5.56 


: -3.65 : 0.35 






9g' ■ 


9f' -- 


= 0.94 : 


2.04 : 6.24 


v 


9% 


■ 9g' 


■ 9s' 


= 0.65 


0.69 : -0.20 



(B2) 



Two-photon decay widths: 

0.897A:e\/ 15^27 



0.095 keV 



0.465 keV T 



0.326 keV 15,^27 = 0.152 keV , 



^N'-*2-y/^ S'-*2-y — 3.05 rG>/_»27/r5/-^.27 — 2.14 



(B3) 
(B4) 
(B5) 



For an even further increase value of the cut-off width A = 2.5 GeV 
KsG = 0.041 GeV-^) we get M^v = 1.287 GeV and Mg = 1.695 GeV. 
Mixing matrix: 

0.78 0.61 0.11 
M= I -0.62 0.74 0.27 
0.08 -0.29 0.95 



0.982 GeV; 



(B6) 



Set of coupling constants: 



/ 


9n' 


9g' 




9f^' 


ss 

9g' 


\ 


9%'^ 


■ 9g' 



• 9s' 

„ss _ 

9s' - 

■ r,99 



= 4.64 : -3.18 : 0.32 \ 
0.79 : 1.70 : 5.27 
= 0.51 : 0.55 : -0.17 / 



(B7) 



Two-photon decay widths: 



rjv-.27 = 0.911 keV , r5_27 = 0.101 keV, 

rAr'-.27 = 0.455 keV rG'-.27 = 0.347 keV rs'-^2'r = 0.165 keV , 

rjv/-^.27/r5/-^27 = 2.99 rG'-*2-y/^S'-^2j — 2.28. 



(B8) 
(B9) 
(BIO) 



As a last point we consider the increase of the effective quark mass parameter from 
0.86 to 1.1 GeV in order to check its influence on the results. For A = 1.5 GeV we get 
fXs = 1-2 GeV and Ksg = 1-00 GeV'^. The masses are Mjv = 1.307 GeV and Ms = 1.695 
GeV. The mixing matrix is 



M 



/0.81 0.57 0.09 \ 
-0.57 0.78 0.25 
0.07 -0.27 0.96 / 



(Bll) 



where again no decisive variation from the previous cases is seen. 
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FIGURES 




q-p/2 



FIG. 1. N- meson mass operator diagram. 




FIG. 2. Glueball mass operator diagam. 
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FIG. 3. N' = /o(1370) meson decaying into two photons through a u-loop. 




FIG. 4. Bubble and tadpole diagrams for the two-photon decay of the N' = /o(1370) meson. 
Both diagrams are suppressed, but conceptually important for conservation of gauge invariance. 
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FIG. 5. Plot of the mixing functions /(p^),r(p^) (eq. (49)) and f*{p^),r*{p^) (eq. (55)) for 
the first propagator choice and with a vertex cut-off A = 1.5 GeV. 
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FIG. 6. Plot of the mixing functions f{p^),r{p^) (eq. (49)) and f*{p^),r*{p^) (eq. (55)) for 
the second propagator choice and with a vertex cut-off A = 1.5 GeV. 
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